Mechano chemical coupling of kinesin studied with the neck linker swing model 
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We have proposed the neck linker swing model to investigate the mechanism of mechanochemical 
coupling of kinesin. The Michaelis-Menten-like curve for velocity vs ATP concentration at different 
loads has been obtained, which is in agreement with experiments. We have predicted that Michaelis 
constant doesn't increase monotonically and an elastic instability will happen with increasing of 
applied force. 

PACS numbers: 87.16.Nn 



I. INTRODUCTION 

Kinesin is a molecular motor that transports organelles 
and membrane-bound vesicles along a microtubule in 
various cells 0, Q- It takes hundreds of 8 nm steps 
(the size of tubulin heterodimers composed of a and 
subunits)HHI3 before detachment and the run length 
is longer than 1 /xm0, [J 0. Q This is why kinesin is 
called a processive motor. 

Conventional kinesin (hereafter called 'kinesin') is a 
dimer consisting of two identical ~ 120kD chains, com- 
monly known as heavy chains. Each heavy chain contains 
a N-terminal globular head domain, a stalk region which 
is responsible for heavy chain dimerization, a tether that 
joins the head and stalk, and a C-terminal fan-shaped tail 
domain which usually binds cargo for transporting in liv- 
ing cell or bead for applying force in single molecule ma- 
nipulated experiments [jllSl]- The head domain is a highly 
conserved region for different members of the kinesin su- 
perfamily and contains two binding sites. One binds to 
microtubule and the other to nucleotides. Two heads al- 
ternately hydrolyze one molecular of ATP for each 8 nm 

ste P [l llllSHHHHElHIil. Sometimes : the 

8-nm step can be resolved into fast and slow substeps, 
each corresponding to a displacement of ~ 4 nm[l4|. The 
tether is a ~ 15-amino-acid segment. It becomes immo- 
bilized and extended towards the forward direction when 
its head binds microtubule and ATP, and reverts to a 
more mobile conformation when phosphate is released 
after ATP hydrolysis p^j. If the tether is replaced by a 
random sequence of amino acids [l6| or if it is cross-linked 
to nucleotides site 01 j the kinesin will lose the capabil- 
ity of stepping. Therefore, the conformational change of 
tether seems to be necessary for kinesin to step|l5|. 

The motility of kinesin can be explained by a "hand 
over hand" model0 El ES EH- The two heads 
alternately repeat single- and double-headed binding 
with microtubule. A simplified binding mode has been 
proposed^. For single- headed binding, the attached 
head either binds ATP or is empty, whilst the detached 
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head binds ADP. For double-headed binding, the forward 
head is empty, whilst the rear head binds either ATP or 
ADP-Pi. 

Kinesin works in a cyclic fashion for several intermedi- 
ate states. The Michaelis-Menten relation for the rate of 
ATP hydrolysis is still a basic law[l|l[lS|2l||l|i||26|. 
The average stepping velocity, however, also has been 
experimentally found obeying the Michaelis-Menten law 
for a range loads (2^, which means the mechanochemical 
coupling is tight, i.e., kinesin hydrolyzes one molecular of 
ATP for each 8 nm step[l2|- With increasing of the ap- 
plied force, the saturating velocity decreases as expected, 
however, Michaelis constant surprisingly increases (2?|. 
The role of force in the reaction kinetics can be used to 
investigate the mechanism of mechanochemical coupling. 

In this paper, we present a neck linker swing model to 
investigate the mechanism of mechanochemical coupling 
of kinesin. We will discuss the effect of applied force 
on the reaction kinetics and the elastic instability of the 
complex composed of neck linker and attached head. 



II. NECK LINKER SWING MODEL 

When a head attaches microtubule, its tether becomes 
a rigid rod that can be bent by an applied force, and we 
call this tether as a neck linker. Otherwise, we call the 
tether as a flexible chain if its head is detachedpjj. 

An idealized scheme of the relations among nu- 
cleotides, microtubule and the tether's conformation is 
shown in Figure Hh[H |H HHU- For attached head, 
if its catalytic cleft is empty, the neck linker tends to be 
perpendicular to microtubule. ATP binding will change 
the catalytic cleft's conformation. The allosteric interac- 
tion between two binding sites in the attached head will 
trigger neck linker to swing to the forward direction 

The neck linker swing model consists of two chemi- 
cal transitions (k\ and k^) and two mechanical substeps 
(&2 and ki) as shown in Figure^). This model is con- 
sistent with the widely accepted model for the kinetic 
mechanism of kinesin |31|. We just keep their main pro- 
cesses and rearrange them into steps in our model. This 
approximation is reasonable for the case that ADP con- 
centration is low. 
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FIG. 1: (a): An idealized scheme of the relations among nu- 
cleotides, microtubule and the tether's conformation(not to 
scale), (b): Neck linker swing model. Two chemical transi- 
tions and two mechanical substeps are coupled alternately to 
each other. T, D and P represents ATP, ADP and P; respec- 
tively. 



the free head and bead will be thrown forward. A 
power stroke occurs. We assume the motion of bead, 
by which optical tweezer can applies force on motor, is 
overdamped, the average velocity is 



<i> = ~«/c) + (/b)-/) 



(2) 



where £ is the viscous coefficient, (/ c ) and (/b) are av- 
erage forces generated from ATP binding and releasing 
of the bending energy respectively, that is, (/ c ) + (/b) = 
(AEi + E^ en( ^)/(^L + £). The power stroke rate then is 



l/C {AE 1+ E hend 



L/2 + e L/2 + i V L/2 + t 



f) (3) 



AEi is the energy that power stroke outputs, and orig- 
inates from the catalytic cleft's conformational change 
induced by ATP binding [Tlj. 



C. chemical transition 2: detached head's ADP 
releasing and attached head's ATP hydrolysis 



In this model, the mechanochemical cycle of kinesin 
includes four steps: ATP binding; power stroke; ADP 
releasing + ATP hydrolysis and recovery stroke. 



A. chemical transition 1: ATP binding 

When a head strongly attaches to the microtubule and 
waits for ATP binding, the neck linker will be bent by 
the applied force / to the backward direction with dis- 
placement £ as shown in Figure^.. The bending energy 
£"b en d ~ and the force / = dE^ end /d£ » kI in the 
case of slight bending, where k is the bending rigidity. 

The attached head will finally bear the bending energy 
and the energy barrier for ATP binding will increase. So 
ATP binding rate can be written as: 



ki(f) = k h e' Eb "" d/kBT [ATP] 



(1) 



where k^ is ATP binding constant without applied force, 
k B is Boltzmann constant, and T is absolute temperature. 
This relationship holds whether [ATP] (ATP concentra- 
tion) is high or low. During attached head is waiting for 
ATP binding, the free head can't reach any binding site 
on the track because the end-to-end distance of the flexi- 
ble chain is not long enough. These sites are separated by 
a fixed distance (L = 8 nm) along the rigid microtubule. 



B. mechanical substep 1: power stroke 

Once the attached head binds an ATP molecule, neck 
linker will move to the new equilibrium position, and 



With starting of the power stroke, the random mo- 
tion of the free head will be biased to the forward direc- 
tion and the flexible chain will be stretched due to neck 
linker's throwing and the interaction from the nearest 
/3-tubulin monomer, which greatly increases the proba- 
bility that the free head reaches the next binding site on 
the track. It first bind weakly, eventually, releases its 
ADP and attaches strongly to the track. The stretched 
chain will shrink and transform into a rigid rod, which 
helps the rear head to split the bound ATP, and one of 
products, phosphate will be released^, |22, H3- This 
reaction weakens the binding of the rear head with the 
microtubule and leads it to detach from the track with 
ADP [l5ll22Tl30l |. We assume the rate of the process from 
the end of power stroke to the detaching of rear head, k$, 
is independent of applied force. 



D. mechanical substep 2: recovery stroke 

The tether linked the forward head now becomes the 
neck linkerplllllil, and will swing to the mechanical 
equilibrium position. A recovery stroke occurs. Similar 
with the power stroke, the rate of this step 
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where AE2 originates from the catalytic cleft's conforma- 
tional change induced by ADP releasing, which is confor- 
mationally the recovery process of ATP binding. 

The cycle is now ready to repeat, the difference is that 
roles of these two partner heads have exchanged. The ki- 
nesin dimer has hydrolyzed one ATP and moved forward 
8 nm with two substeps [H [H IHH- Although there 
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are more than two chemical transitions as above, all other 
chemical processes are rapid rate transitions (random- 
ness shows that there are only two to three rate-limiting 
transitions |2H l32l l33| ) . We lump all other chemical tran- 
sitions into &3. 



III. FITTING AND RESULTS 

From the time series shown in Figure^), the average 
time to complete a single enzymatic cycle, (r), can be 
computed conveniently |33J, |34j. The average velocity of 
kinesin moving along microtubule is 
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where v max and Km are saturating velocity and Michaelis 
constant respectively. Obviously, the average velocity 
obeys the Michaelis-Menten law as observed [27|. 

In our model, AE\ and AE2 are independent of ATP 
concentration. We define the efficiency of mechanochem- 
ical coupling of kinesin as r\ = (AEi + AE 2 )/Afj, and as- 
sume it is a constant, where A/i is the free energy excess 
in one molecule ATP hydrolysis. It must be noted that 77 
defined here is different from the efficiency of motor. We 
assume kinesin has adjusted AEi and AE2 to achieve 
optimal kinetic velocity in evolution. From Eq.@ and 
Eq.gJ), we yield AEi = AE 2 = AE (< 5 A/i), which 
is reasonable with the fact that the catalytic cleft's con- 
formational change induced by ADP releasing is confor- 
mationally the recovery process of that induced by ATP 
binding. For simplicity, we introduce the rate of power 
(or recovery) stroke without load, ko = AAE /((X 2 ). 

According to Eqs.® and (@J, the power stroke rate k 2 
decreases with increasing of applied force /, whilst the 
recovery stroke rate k^ increases as shown in FigunJ2Kd). 
At a critical applied force, k 2 becomes zero. We call this 
critical force as stall force. It can be derived as 



/stall 



AE + i^bend (/stall 

L/2 + £(/ at all) 



(8) 



If applied force is larger than the stall force, the force in- 
duced by conformation change can't overcome it to push 
the kinesin forward along the microtubule. 

For the case of slight bending of neck linker, Eq.JSJ) 
can be rewritten as / Btall « 2AE/[L{\ + f staU /{nL))} < 
A/J./L ~ 10 pN with A/i w 80 pN-nm under physiolog- 
ical conditions 1351 . which is in agreement with the ex- 
perimental data 27] . It must be noted that stall force in 



TABLE I: The fitted parameters for experimental Umax and 
Km versus f|27| with k B T r — 4.1pN-nm 
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6.50 1.88 1.33 
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1.45X10" 3 



our model is independent of ATP concentration. We can 
rewrite Eq.© and Eq.JTJl as 
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(10) 



and use them to fit the measured data of w m ax and K m at 
different loads^, where f A = //[/ stall (l + f stall / (kL))}, 
and / = f /(kL). The fitted values of / ata ii, k, k^, ks and 
( are listed in tabled and the two fitted curves, w m ax ~ / 
and Km ~ /, are shown in Figure [2K a). 

The fitted stall force / Btall = 6.5 pN is consistent with 
measured value|27|. The total energy outputted in one 
cycle, 2A_E w 74 pN-nm, is less than A/x, which demon- 
strated again that the fitted values are reasonable. The 
efficiency of mechanochemical coupling of kinesin, 77, is 
about 93%. However, the efficiency of motor, / atElll L/A/i, 
is about 65%, which is in agreement with experiments. 

It is the motion of a silica bead that is measured in 
experiment [23]. The radius of the bead, R, is 0.25 /zm. 
The viscousity of water at room temperature, fj w , is 
about 0.9 x 10~ 9 pl\l-s/ nm . The Stokes drag coefficient 
of the bead can be estimated by 6irRf) w , and is about 
4.24x 10 -3 plM-s/nm. The fitted viscous coefficient of bead 
in our model, thus, is comparable with the Stokes drag 
coefficient for a sphere in 0.5 /tm diameter at room tem- 
perature. 

The widely accepted model for the kinetic mechanism 
of kinesinj^l has proposed that the ATP binding rate 
is about 2 fiM~ 1 -s~ 1 and the rate from free head's ADP 
releasing to attached head's detaching can be estimated 
about 70 s _1 . In our model, ATP dissociating rate k_\ 
isn't taken into account, this is why k^ is slightly less than 
the proposed value. k% also approaches the estimated 
value. 



IV. DISCUSSION 

A. no longer increases monotonically with load 

Although the recovery stroke will go faster and faster 
with increasing of applied force, the power stroke will 
spend more and more time and finally stalls at the stall 
force, which leads to the concave down of saturating 
velocity. With low applied force, Michaelis constant 
will nearly increase exponentially because e Ebcnd ^ kBT in- 
creases faster than v mSLX decreases. With high force, 
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FIG. 2: (a): Fitting to experimental Umax and Kyi vs f |27j|- 
The fitted parameters are listed in table U] (b): Dependence 
of velocity on load at saturated and limiting ATP concentra- 
tions. The curves are directly from Eq.JSJ analytically with 
parameters shown in tabled and the experimental data come 
from Ref 27]. (c): Applied force vs extension. The relation 
/ = k£ with k shown in table dis in good agreement with the 
experimental data 22] at low force, (d): Dependence of rates 
on applied force at 0.1 mM ATP. 



power stroke becomes the rate-limiting step for the sin- 
gle enzymatic cycle, so, Ky[ will fall fast as v mSuX shown 
in Figure^ a). Michaelis constant no longer increases 
monotonically with load. The big error bar at high force 
implies that this prediction is possible. 



with increasing of applied force, and its force-extension 
relation can be fitted by a simply elastic model again. 



C. velocity versus applied force at saturating and 
limiting ATP concentration 

With these reasonable fitted parameters, we can use 
Eq. (JSJ directly to compute the dependence of the velocity 
on load at saturating and limiting ATP concentrations. 
As shown in Figure Eljb) , the theoretical velocity-force 
relation is in good agreement to experimental data|27|. 
It is clear that there are three distinctive regimes of ap- 
plied force: (1) If the load is low, chemical transition 2, 
the process from free head's ADP releasing to attached 
head's ATP hydrolysis, k$, is the rate-limiting transition 
at saturating [ATP] and v w Lk^ rj 900 nm/s, which is 
consistent with what is known about the biochemistry 
of kinesin[3^.l3e|. while ATP binding is the rate-limiting 
transition at low [ATP] and v w Lk% » 50 nm/s with 
5 /iM. (2) If load is high, kinesin will stall at the same 
force as discussed in section IIIII whether ATP concen- 
tration is high or low. (3) If applied force is moderate, 
the velocity-force curves display different shapes at low 
and saturating ATP concentration respectively. At very 
low ATP concentration, the velocity decreases exponen- 
tially as k\ and looks like linear with load because ATP 
binding is the rate- limiting step. At saturating [ATP], 
the velocity is concave down with load as discussed in 
section llVAl 



D. two mechanical substeps 



B. neck linker's bending 

Actually, the neck linker's bending has been uncon- 
sciously revealed in experiment 22]. With an applied 
force moving toward the end of microtubule, a simply 
elastic model can fit all extension of kinesin-microtubule 
complex for double-headed attaching. But the force- 
extension relation in Figurt^c) can't be fitted by a same 
simply elastic model for all extension for single-headed 
attaching. With our neck linker swing model, we can 
image there are three distinctive regions of force. If the 
force is low, the "Extension" is mainly contributed by 
neck linker bending. We compare the relation f = k£ 
to the experimental force-extension data 22] in Figure 
121(c). It is very surprising that our model is in agree- 
ment with experiment very well. With the increa sing 
of applied force, a phenomenon of elastic instability |3a] 
will inevitably happen to the complex composed of neck 
linker and attached head. This is why the extension 
rapidly increases while the "applied force" falls fast as 
the measured data in Figure[2(c) . After reaching the new 
mechanical equilibrium, the complex will further extend 



The two mechanical substeps in this model are con- 
tributed by the two heads respectively. The recovery 
stroke is always a rapid rate step which corresponds to 
the observed fast substep[3- If a moderate load such as 
1.5 pN acts on the bead and ATP concentration is main- 
tained at 0.1 mM, the two chemical transitions have the 
same rate as shown in Figuri^d). The time spent in 
ATP binding equals that spent in chemical transition 2. 
Theoretically, 4.8 nm slow substep and 3.2 nm fast sub- 
step can be detected directly by single molecular manip- 
ulated techniques such as optical tweezers. The different 
load-dependence of these substeps' rate may be revealed 
in the future experiment. 



V. CONCLUSION 

We proposed the neck linker swing model which divides 
the single enzymatic cycle into two chemical transitions 
and two mechanical substeps. Each chemical transition 
will induce the conformational change in the catalytic 
cleft and generate a corresponded mechanical stroke. The 
model can be used to explain the observed substeps [Tij]. 
The different load-dependence of these two strokes' rate 
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may be revealed in the future experiment. We have in- 
vestigated the mechanism of mechanochemical coupling 
of kinesin by the influence of applied force on the bend- 
ing of neck linker. When attached head is waiting for 
ATP binding, the neck linker is bent by the applied force. 
The attached head bears the neck linker's bending energy 
and the energy barrier for ATP binding increases. This is 
why Michaelis constant increases with applied force. Our 
theoretical analysis of average velocity of motor in Eq. JSJ) 
also obeys Michaelis-Menten law and has been used to fit 
the observed saturating velocity and Michaelis constant 
at different loads j^]. The fitted values of chemical re- 
action rates are in agreement with those in the widely 
accepted model and the fitted viscous coefficient of 
bead is also comparable with the Stokes drag coefficient. 
The stall force is independent of ATP concentration and 
its fitted value is consistent with the observed data in 
experiment The fitted bending rigidity of neck linker 
can be used to explain the relation of force-extension in 
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